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Professor H. G. Seeley's objection to the importation of the Amphi- 
bian plan in explanation of a part of the Plesiosaurian shoulder girdle, 
of which other parts have been explained by reference to the Chelonian 
plan of construction, has not, I venture to think, great weight, since, 
of the early Reptilia, from the time when their remains first began 
to be studied, it has been a frequent remark that their skeletons 
comprise structural arrangements which, in existing animals, are now 
found separately. Moreover, it is thought by some of the ablest 
comparative anatomists that the Chelonian skeleton shows closer 
approach to the Amphibian than is to be found elsewhere. 
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(Abstract.) 

1. The object of the present paper is to show how to determine 
expressions for the current in circuits having iron cores, similar to 
the well-known equations 

*' = JL(l_< r B*/i') 

R v J 

for circuits without iron, which will enable the current curves to be 
pre-determined by calculation and plotted independently of experi- 
ment. 

In circuits with iron cores the value of — - occurring in the original 

dt & & 

differential equations 

coxy -rj » /i \ 
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dt 



J.J bill l/f —♦ ZZZ JLvx>* » ..4. ..««<,. .6... i & ] 

dt 

continually alters as i changes. If we could obtain an expression for 

dB 

~~~ in terms of i, the substitution of this expression in (1) and (2) 

should lead us to the required result. But, though such an expres- 
sion can be found, its substitution will generally lead to differential 
equations which cannot be solved by known methods. 

2, In the case represented by (1), where the applied E.M.F. is 
constant* we can determine by Lagrange's formula of interpolation 
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tlie equation to the (B, H) curve of the particular core under 
consideration. This will he of the form 

B = a -f ajl -f a 2 H 2 -f- . . . . +a»H w (3), 

where n is one less than the number of observed simultaneous values 
of B and H from which the equation is calculated ; whence 

dB T cli , T2 .di , . T»«»-i^ /,<\ 

— -— = «iJu— -f 2a 2 m — j- .... -\-na n hH >l ■ — • (4), 

dt dt dt dt 

and substituting in the equation 

we get 

a 1 L-f-2a 2 L 2 ^H-3a 3 L 3 r+ . . . . -j~ na n LH u ~" 1 



E-Bi 



cZ* = cZ£ . . . . (6), 



which is easily integrable ; and integrating between the proper 
limits we get after reduction 

E / l E E 2 E"- 1 



, l E E 2 E H ~ 2 

/ l E E 2 E w ~ 3 \ 

™i^(^3a 3 L J -J-^4a 4 LH-T^5a 5 L 4- • • • • + g-^ wa>,L w ] 

/ 1 E E 2 E w - 4 

— £ z 3 -^ 4 %L 4 + =j5 5 a 5 L 5 + = 6 a 6 L 6 + ....+ — -- wa^L' 

\K K~ li° li M ~ 3 

The corresponding equation when the E.M.F. is removed and the 
current is dying away is 

, / iL 1 E 



R °Il 



& 



+ R 2/ ~ R+^ L ^ + ^ L + .... + w/jl L 

/ 

-■^(a/ 3 L a t + i3/ 8 L 8 ^ + i4/ 4 LV+ . . . . +^j«/«L w * M " 1 ) • • (8), 

the fu f 2, f 3 9 &c., being the coefficients of the powers of II in tbe 
equation to the descending (B, H) curve, which is, of course, different 
to the ascending curve. 
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3. This method is not applicable to the case in which the impressed 
B.M.F. is sinusoidal, on account of difficulties of integration. But 
both cases can be treated in another way : — Take a series of points 
on the (B, H.) curve of the iron core, such that the chords joining them 
practically coincide with the curve itself. Let B K , H, c and B K +i, H K +! 
be the coordinates of two consecutive points. The equation to the 
curve between these points is approximately 

B = )w K +! EL -j- constant (9), 

where m K+1 = ~- _, 

and therefore between these limits 

(a/JL* tt Ctl/ /T a\ 

dt " dt 

During the time that the current rises from i K to i K +i 9 and B and 
H rise from B* and H* to B* + i and H^ + i, and t rises from t K to t K + i t 
we have 

di 



and therefore 



4 - J. I m K + l^ , E — Ri K HON 

r K+1 — c K -| - — • log ,., yy . — • • • k-l^Ji 



which is true to a very close approximation for any simultaneous 
values of t and i between the above limits. From this equation, 
since t and i are both zero, we can determine in succession the times 

at which the current has the known values 0, — -, — ?, . e . # &c, 

L L 

using that value of m which applies to that particular value of H 

under consideration. In this way the current curve can be plotted. 

On making E = in the original differential equation, and 

observing the proper limits, we get 

/ _/ i ™* +lL w Hn a*} 

bn+i — hiiT -q 1U & vr \ x '- J 

as the equation to the curve representing the dying away of the 
current when the E.M.F. is withdrawn; m m m n+ i being determined 
from the descending (B, H) curve. 

Fig. 1 and Table I give the results of calculation for a circuit with 
the following constants : — Resistance, 1 ohm ; E.M.F. ? 0'4315 volt ; 
self-induction (without iron core), O^OOO^ henry. 
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Table J, 



Ascending. 


Descending. 


Current. 


Time. 


Current. 


Time. 


-oioo 

-0200 
0-0270 
-0460 
0*0600 
-0825 
0-1075 
0-1400 
0-1905 
-2635 
0-3545 
0-8750 
-3875 
-4000 
-4125 
'4315 


-00023 
-00095 
-00228 
0-01025 
-01649 
-02269 
-02696 
-03100 
0-03456 
-03768 
-04048 
'04150 
-04205 
-04278 
0-04367 
oc 


-43150 
-26350 
0-19050 
0-14050 
0-08250 
-04600 
-02700 
-01775 
-00925 
0-00475 
-00125 

o -ooooo 

1 


'000000 
-000235 
-000768 
-001460 
-003036 
-007670 
-020009 
039880 
-067370 
-101070 
0-185770 
a 




4. When the impressed E.M.F. is sinusoidal, we substitute for 
dB/dt in the equation 

dB 



E sinpt- 



dt 



-Lvi/ ••••••••«•««»##• IJL'^ri) 



having determined the various values of dB fdt, as in the foregoing. 
As by the present method the value of m changes abruptly from 
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m K to m K +u w ^ must employ the general solution of (14), wliicli for 
the interval t KJ t K + x is 



y(R 2 + m K+i yL 2 ) 

in order that the current at the commencement of the interval 
t K , t K +i may have the same value which it had at the end of the 
interval t KM \t K ^ The complementary function 

enables us to ensure this condition ; for, by taking the constant 
A K +i of such a value that equation (15) is satisfied when i = i K and 
t = t K , there is no abrupt change in the current. The complementary 
function, in fact, represents the gradual dying away of whatever 
excess or defect of current there would be in the circuit when m 
changes. 

Equation (15) is true for all values of i between i K and i K +\ ; and, 
therefore, enables us to find the time £ K+1 at which the current 
attains the known value H K + 1 /L. 

By changing k into k + 1 we obtain similarly the time t K+2 at 
which the current has the value H k+ 2/Ij, and so on. 

Thus the determination of t K+1 is made to depend upon t K9 and in 
order to make a start we must assume that the value of i is known 
for some definite value of t. When the number of alternations 
per second is not great, it is not of much consequence what assump- 
tion, within reason, is made, as, though the calculated curves will 
vary with the assumption made, they will all eventually merge 
into the true periodic current curve at some point which will be 
exhibited when the first evanescence of Ae~ R '/ mL takes place. 

As this complementary function is a continually decreasing quan- 
tity, it becomes negligible when it is allowed time enough. This 
opportunity is afforded when the straighter portions of the (B, H) 
curve are reached, and where the points on the curve can be taken 
further apart. 

When, however, the period of alternation is short in comparison 
with the time-constant of the circuit, the evanescence of Ae"~ E '/ mL 
< does not so readily take place ; and it will generally take several 
: cycles before the current shakes down into its truly periodic forni. 
The preliminary assumption ought 'therefore in such cases to be 
made with care if it is desired to avoid the labour of calculating the 
first cycles. But, if the periodicity is quick enough, Ae~ R ' /wlL can be 
taken as a constant, at any rate during the shorter intervals. 

Fig. 2 gives the plotted curve calculated for a circuit consisting of 
500 turns surrounding an anchor ring, having a coefficient of self- 
induction (without the core) of 0'0004j henry, and a resistance of 

s 2 
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1 ohm. The impressed E.M.F. is 12*5 volts. The periodicity is 
slow, being 6 cycles per second, and the true periodic curve is seen 
to appear before the end of the first half period. Curve A is calcu- 
lated with the preliminary supposition that i = when t = 0, and 
Carve B on the assumption that i — when ft = 0, neither of which 
is strictly correct. 

Fig. 3 gives the (B, H) curve, and Table III the observations 
upon which it is based, of the iron core under consideration. 

5. The method of this paper can be applied to the pre -determination 
of the curves of the primary and secondary currents in transformers 
and of the curve of magnetic induction with regard to time. 

[For, when there is supposed to be no magnetic leakage in the core, 
the expressions for the primary and secondary currents are- — 

Ey(VS 2 +|) 2 mVP) . 

x = ==—■ sin( pt— 0)-]-Ae $* . ... (lb), 

y = ' cos (j)£— 0)-f Be £' (17), 

, , . nfpmhS 2 

where tan 6 = 



% 4 RSH%ym 2 L 2 (% 2 R4-^i 2 S) ' 
, wmL(^ 2 2 R + ni 2 S) 

D = v/{% 4 R 2 S 3 +p 2 m 2 L 2 0? 2 2 R + % 2 S) 2 } 

_ _ n?R&_ 

* ~" mL0i 2 2 R-f*% 2 S) ' 

R being resistance of primary, S that of secondary, n x the nnmber of 
primary turns, n 2 the number of secondary turns, L the self-induction 
of the primary. 

From (16) and (17) H, the total magnetic force on the core (being 



*^±M L \ reduces to 
% 2 / 



?!i^E sin (irf-0) +Ce-*' (18), 

1) 



where v - %2 



tij.A4-%B r 

O = ~7~o. U ' 



Now equation (18) treated in the same way as (15) gives the 
simultaneous values of H and U that is, of m and t. These latter 
substituted in (16) and (17) give the simultaneous values of x and t 
and y and t. 
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If there is magnetic leakage it is necessary to discriminate between 
the H of the primary core and the H of the secondary core ; these 

are respectively-— 



Hi 



and H 



E y{L 2 S 2 +i)"m 2 (LrT-M 2 ) 2 } sin (pt-e) 

n x * v / [{RS- j p 2 m 8 (LN--M 2 )} 2 +jp 2 m 2 (NR-|-LS) 2 ] 

+ Cie"^'+ C 2 e- f Jj ...*.. (19), 

E MS sin (pi—ty 



~ wa v /[{RS-j>W(LN~M 3 )} a +jp ? m s (I i rR + LS) 2 ] 

+ Die"^H-'D 2 e"«' 9 < ...... (20), 

derived respectively from H = "— (La? + My) and H 2 = — (Ma? ■+• ~Ny) 9 

E */(S 2 -f p 3 m 2 N 2 ) sin (pi— ?/) , . . . , ,„... 

where a? == --^^Zlii^^ vi 7 '_Z + A^-^'-f A 2 e-* 9 ' 8 * , . (21), 

and y = — _£^L_?~l£JIlili ^-Bje-^-f B 2 e~-'M ......... (22). 

In these equations 1ST is the self- induction of the secondary, M is 
the mutual induction, F is the radicle in the denominator of (19) and 

(20), wliiie tan -« = pm — — — ^ ~-^— , 

1 /J ; ^ R(S 3 +p 5i m 2 lN 3 )-f_p 2 m 2 M 2 S 9 



tan Y r = JP"^ 



RS-jp^LN-M 2 ) 

, L 2 S 3 -fRSM 2 -f-^m 2 (LF-M 2 ) 2 
tan g = »?n, — _— ±- ^ — . — dL , 

^ R{LS a +p 2 m*N(LN -M 2 ; } ? 

___ ^R;+L S^-y/{(NR-LS) 2 + 4RSM 3 } 

qi " 2m(TlF~M 2 ) 

= rTR + LS + v / {fflR-LS) 2 -f4RSM 2 } 
% 2m(LN-M SJ ) 
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